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Sur la stabilite et I'approximation de flots de distribution de 
processus de branchements, avec applications au filtrage non 

lineaire multicibles 



Resume : Nous analysons les proprietes de stabilite exponentielle d'une classe d'equations 
a valeurs mesures que Ton rencontre dans des problemes de filtrage non-lineaire multicibles. 
Nous demontrons ensuite les proprietes de convergence uniforme, par rapport a I'horizon 
temporel considere, d'une famille assez generale d'algorithmes de filtrage multicibles. Cette 
analyse s'applique notamment aux methodes de type Monte Carlo sequentielles et aux algo- 
rithmes particulaires fondes sur revolution de systemes de particules en interaction de type 
champ moyen. Nous illustrons ces resultats dans le cadre des filtres de Bernoulli et les filtres 
PHD (Probability Hypothesis Density). Ces resultats scmblcnt etre les premiers de ce type 
pour ces classes de modeles de filtrage stochastiquc multicibles. 

Mots-cles : Processus a valeurs mesures, filtrage non-lineaire multicibles, filtre de Bernoulli, 
filtre PHD (Probability hypothesis density), filtres particulaires, systemes de particules en 
interaction, techniques dc champ moyen, proprietes de concentration exponentielle, inegalites 
de contraction fonctionnelles. 
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1 Introduction 

Let {En)n>o be a sequence of measurable spaces equipped with the cr- fields {£n)n>o, and for 
each with n > 0, denote M.{En), Ai+{En) and V{En) the set of all finite signed measures, 
the subset of positive measures and the subset of probability measures, respectively, over the 
space En- The aim of this work is to present a stochastic interacting particle interpretation 
for numerical solutions of the general measure- valued dynamical systems 7^ G A^+(£'„) 
defined by the following non-linear equation 

7„((ix„) = (7„_iQn,^„_J (dxn) := / 7n-i(d2;„_i)(5„,^,^_-^ (x„_i, dx„) (1.1) 

with initial measure 70 G A4^{Eq), and positive and bounded integral operators Qn,^ from 
into En, indexed by the time parameter n > 1 and the set of measures 7 G 

This class of measure- valued equations arises in a natural way in the analysis of the first 
moments evolution of nonlinear branching processes, as well as in signal processing and more 
particularly in multiple targets tracking models. A pair of filtering models is discussed in 
some details in section [1.1. II and in section [1.1. 21 In the context of multiple targets tracking 
problems these measure-valued equations represents the first-order statistical moments of 
the conditional distributions of the target occupation measures given observation random 
measures obscured by clutter, detection uncertainty and data association uncertainty. 

As in most of the filtering problems encountered in practice, the initial distribution of 
the targets is usually unknown. It is therefore essential to check wether or not the filtering 
equation "forgets" any erroneous initial distribution. For a thorough discussion on the 
stability properties of traditional nonlinear filtering problems with a detailed overview of 
theoretical developments on this subject, we refer to the book [5] and to the more recent 
article by M. L. Kleptsyna and A. Y. Veretennikov [12J. Besides the fact that significant 
progress has been made in the recent years in the rigorous derivation of multiple target 
tracking nonlinear equations (see for instance [H [21 |2U [E]), up to our knowledge the 
stability and the robustness properties of these measure-valued models have never been 
addressed so far in the literature on the subject. One aim of this paper is to study one such 
important property: the exponential stability properties of multiple target filtering models. 
We present an original and general perturbation type technique combining the continuity 
property and the stability analysis of nonlinear semigroups of the form (jl.ip . A more 
thorough presentation of these results is provided in section 11.21 dedicated to the statement 
of the main results of the present article. The detailed presentation of this perturbation 
technique can be found in section [3l 

On the other hand, while the integral equation (jl.ip appears to be simple at first glance, 
numerical solutions are computationally intensive, often requiring integrations in high di- 
mensional spaces. One natural way to solve the non-linear integral equation (jl.ip is to use 
find a judicious probabilistic interpretation of the normalized distributions fiow given below 

r]n{dxn) := 7n(c^a:„)/7„(l) 

To describe with some conciseness these stochastic models, it is important to observe that 
the pair process (7„(l),7/„) S (M+ x V{En)) satisfies an evolution equation of the following 
form 

(7n(l),?/n) = r„(7„_i(l),r/„_i) (1.2) 
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Let the mappings : M+ x V{En) M+ and : x V{En) V{En), denote the first 
and the second components of r„ respectively. By construction, we notice that the total 
mass process can be computed using the recursive formula 

7n+i(l) = 7n(G„,^„) = ??„(G„,^,J 7„(1) with G„,^„ := (1.3) 

Suppose that we are given an approximation (7^ (1), 7/^) of the pair (7^(1), ?]„) at some 
time horizon n, where N stands for some precision parameter; that is (7^ (1), 77^) converges 
(in some sense) to (7,1(1), 77^), as — >• 00. Then, the A''-approximation of the measure 7„ is 
given by 7^ = 7^(1) x rj^ . The central idea behind any approximation model is to ensure 
that the total mass process at time {n + 1) defined by 

l^+i{^)=ri'^.{Gna-)ln{^) (1-4) 

can be "easily" computed in terms of the iV-approximation measures 7^. Assuming that the 
initial mass 7o(l) = 7(^(1) is known, the next step is to find some strategy to approximate 
the quantities r,^_^]^(7^ (1), ry^ ) by some A/"-approximation measures r/^^, and to set 7^1 = 
7^+1(1) xr/^+i. 

The local fluctuations of rj^ around the measures r^(7^]^(l), r/^_^) is defined in terms 
of a collection of random fields : 

W!! := [n^ - r^(7.^_i(l),ry^_i)] ^ ry^ = {l^.A^U^-i) + ^ (1.5) 

which satisfies for any r > 1 and any test function / with uniform norm ||/|| < 1, 

IE(T4^„^(/)|-F„^_i)=0 and E (|M^„^(/)f | ' < a, (1.6) 

where J^^-i — ^ {''Ip ^ < p < nj is the ci-field generated by the random measures rj^ , < 
p < n, while h and ar are universal constants whose values do not depend on the precision 
parameter N . The stochastic analysis of the resulting particle approximation model relies 
on the analysis of the propagation of the local sampling errors defined in p.5|) . The main 
objective is to control, at any time horizon n, the fluctuations of the random measures 
{in ^Vn) around their limiting values (7n,??n) defined by the following random fields: 

K''^ := ^ [7^ - 7n] with K"'^ := ^ K - ^n] ■ (1.7) 

The construction of the A^-approximation measures r]^ is far from being unique. In the 
present article, we devise three different classes of stochastic particle approximation models. 
These stochastic algorithms are discussed in section [H The first one is a mean field particle 
interpretation of the flow of probability measures r/„, and it is presented in section I^Tl The 
second model is an interacting particle association model while the third one is a combination 
of these two approximation algorithms. These pair of approximation models are respectively 
discussed in section 14.21 and in section 14.31 In the context of multi-target tracking models, 
the flrst two approximation models are closely related to the the sequential Monte Carlo 
technique presented in the series of articles [201 IHl ESI IMl 123 [2H], and respectively, the 
Gaussian mixture Probability Hypothesis Density filter discussed in the article by B.-N. Vo, 
and W.-K. Ma |221I23]. and the the Rao-Blackwellized Particle multi-target filters presented 
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by S. Sarkka, A. Vehtari, and J. Lampinen in [n\ I18j. These modern stochastic algorithms 
are rather simple to implement and computationally tractable, and they exhibit excellent 
performance. 

Nevertheless, despite advances in recent years O [TTl I21j . these Monte Carlo particle 
type multi-target filters remain poorly understood theoretically. One aim of this article is 
to present a novel class of stochastic algorithms with a refined analysis including uniform 
convergence results w.r.t. the time parameter. We also illustrate these results in the context 
of multi-target tracking models, yielding what seems to be the first uniform results of this 
type in this subject. 

The rest of the article is organized as follows: In section 11.11 we illustrate the ab- 
stract measure- valued equations (jl.ip with two recent multi-target filters models, namely 
the Bernoulli filter and the Probability Hypothesis Density filter (abbreviate PHD filter). 
Section 11.21 is devoted to the statement of our main results. In section [21 we describe the 
semigroups and the continuity properties of the nonlinear equation 11.11 We show that this 
semigroup analysis can be applied to analyse the convergence of the Bernoulli and the PHD 
approximation filters. Section [3] is devoted to the stability properties of nonlinear measure- 
valued processes of the form ()1.2p . We present a perturbation technique and a series of 
functional contraction inequalities. In the next three sections, we illustrate these results in 
the context of Feynman-Kac models, as well as Bernoulli and PHD models. Section|3]is con- 
cerned with the detailed presentation and the convergence analysis of three different classes 
of particle type approximation models, including mean field type particle approximations 
and particle association stochastic algorithms. Finally, the appendix of the article contains 
most of technical proofs in the text. 

1.1 Measure- valued systems in Multi-target tracking 

The measure-valued process given by (jl.ip is a generalisation of Feynman-Kac measures. 
Its continuous time version naturally arise in the modeling and analysis of the first moments 
of spatial branching process O [8j . 

Our major motivation for studying this class of measure- valued system stems from ad- 
vanced signal processing, more specifically, multiple target tracking. Driven primarily in 
the early 1970's by aerospace applications such as radar, sonar, guidance, navigation, and 
air traffic control, today multi-target filtering has found applications in many diverse disci- 
plines, see for example the texts [T], f2] [T5] and references therein. These nonlinear filtering 
problems deal with jointly estimating the number and states of several interacting targets 
given a sequence of partial observations corrupted by noise, false measurements as well as 
miss-detection. This rapidly developing subject is, arguably, one of the most interesting con- 
tact points between the theory of spatial branching processes, mean field particle systems 
and advanced signal processing. 

The first connections between stochastic branching processes and multi-target tracking 
seem to go back to the article by S. Mori, et. al. [TB] published in 1986. However it was 
Mahler's systematic treatment of multi-sensor multi-target filtering using random finite sets 
theory \10\ [9l [T3l [T3] that lead to the development novel multi-target filters and sparked 
world wide interests. To motivate the article, we briefly outline two recent multi-target 
filters that do not fit the standard Feynman-Kacs framework, but fall under the umbrella 
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of the measure- valued equation (jl.ip . The first is the BernoulH filter for joint detection and 
tracking of a single target while the second is the Probability Hypothesis Density filter. 



1.1.1 Bernoulli filtering 

A basic problem in target tracking is that the target of interest may not always be present 
and exact knowledge of target existence/presence cannot be determined from observations 
due to clutter and detection uncertainty [X5j. The Bernoulli filter is a generalisation of 
the standard Bayes filter, which accommodates presence and absence of the target [23]. In 
a Bernoulli model, the birth of the target at time n + 1 is modelled by a measure ^n+i 
on En+i- The target enters the scene with a probability < 1 and its state is 

distributed according to the normalised measure /x„+i/;U„+i(l). At time n, a target Xn has 
a probability Sni^n) of surviving to the next time and evolve to a new state according to a 
given elementary Markov transition M„_|_i from £"„ into E^+i- At time n + 1, the target (if it 
exists) generates with probability (i„+i(X„+i) an observation l^+i on some auxiliary state 
space, say E^_^^ with likelihood function y). This so-called Bernoulli observation 

point process is superimposed with an additional and independent Poisson point process with 
intensity function /i„ > to form the occupation (or counting) measure observation process 

In its original form, the Bernoulli filter jointly propagates the probability existence of the 
target and the distribution of the target state [23] • Combining the probability of existence 
and the state distribution into a single measure, it can be shown that the Bernoulli filter 
satisfies the integral equation (jl.ip . with the probability of existence of the target given by 
the mass 7n(l) and the distribution of the target state given by the normalised measure 
Vn = 7n/7n(l)- The integral operator for the Bernoulli filter takes the following form 

^ / J N SniXn)gniXn)Mn+i{Xn,dXn+l) + {'Jniiy^ - l)lJ'n+l{dXn+l) 

(1 -7n(l)) +7n(5n) 

where gn is a likelihood function given by 

9niXn) ■■ = i'i- - d„{Xn)) + dniXn)ynilniXnr)/hn) (1-9) 



1.1.2 PHD filtering 

A more challenging problem arises when the number of targets varies randomly in time, 
obscured by clutter, detection uncertainty and data association uncertainty. Suppose that 
at a given time n there are targets (X* ) 

i<i<N^ each taking values in some measurable 
state space En- A target X^, at time n, survives to the next time step with probability 
Sn{Xn) and evolves to a new state according to a given elementary Markov transition M^_^-^ 
from En into En+i- In addition X^ can spawn new targets at the next time, usually 
modelled by a spatial Poisson process with intensity measure Bn+i{Xn, •) on At the 

same time, an independent collection of new targets is added to the current configuration. 
This additional and spontaneous branching process is often modeled by a spatial Poisson 
process with a prescribed intensity measure Hn+i on En+i- Each target X^_^_^ generates 
with probability (i„+i(X*_,_^) an observation Y^^-^ on some auxiliary state space, say -E^+i, 
with probability density function gn+i{X^_^i,y). In addition to this partial observation 
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point process we also observe an additional and independent Poisson point process with 
intensity function Multi-target tracking concerns the estimation of the random measures 
Xn+i = J2i<i<N^ ^x^^ given the observation occupation measures = J2i<i<N"^ ^Y^- 

The multi-target tracking problem is computationally intractable in general and the 
Probability Hypothesis Density PHD (filter), is an approximation that propagates the first- 
order statistical moment, or intensity, of the multi-target state forward in time |14] . The 
PHD filter satisfies the integral equation (jl.ip . with the integral operator given below 

+ 7„(1) I^Ln+lidXn+l) (1-10) 

where Mn+i is a Markov kernel defined by 

)M^_^^-^{Xn, dXn+l) + Bn+l{Xn 1 + 1 ) /-I -I -1 \ 

Mn+l{Xn,dXn+l) ■= 7 \ , h ( ^ ^ 

with the branching rate hn{xn) = The likelihood function (7n,7„ is given by 

5'n,7„ := Tn X 5„,^„ with r„ := {sn + bn) (1.12) 

and 

yn{dy) 7—7 — X V^^' t vT (1-13) 

K{yn) +ln[dngn{-,yn)) 

Since its inception by Mahler [13] in 2003, the PHD filter has attracted substantial interest 
to date. The development of numerical solutions for the PHD filter [21], [23] have opened 
the door to numerous novel extensions and applications. More details on the derivation of 
the PHD filter using random finite sets, Poisson techniques or random measures theoretic 
approaches can be found in the series of articles [H dU [TO]. 



1.2 Statement of the main results 

To describe with some conciseness the main result of this article, we need to introduce some 
notation. We let Osci (£"„), be the set of iSn-measurable functions / on En with oscillations 
osc(/) = sup^ ^/ |/(a;) — /(a;')| < 1. We denote by /i(/) = / fi{dx) f{x) the Lebesgue 
integral of / w.r.t. some measure /i G M.{En), and we let \\^ — v\\tY be the total variation 
distance between two probability measures u and /x on En- 

We assume that the following pair of regularity conditions are satisfied. 

{Hi) : There exists a series of compact sets In C (0, oo) such that the initial mass value 
7o(l) G Iq, and for any m £ In V ^ T^iEn), we have 

6-^n{fn) < rj{Gn,mr}) < G+,n{'ni) for some pair of positive functions 



The main implication of condition {Hi) comes from the fact that the total mass processes 
7n(l) and their A^- approximation models 7^ (1) are finite and they evolves at every time n 
in a series of compact sets 

/„, C [mn,mn] C (0,00) 
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with the sequence of parameters m„ defined by the recursive equations w.„+i = m„ 6'_^„(m„ ) 
and m^^^ = m+^_(_^„(m+), with the initial conditions niQ = tUq = 7o(l)- 

{H2) : For any n > 1, f £ Osci (£'„), and any {m,r]), {m',rj') G (/„ x V{En)), the one 
step mappings r„ = (r^,r^) defined in satisfy the following Lipschitz type inequalities: 

\rl^{m,r]) - r^(m',77') 

\[rl{m,r,)-rUm',7i')] (/) 

for some finite constants c{n) < 00, and some collection of bounded measures S,^ and 
^l,{m',r,')(f^ •) t3{E„) such that 

j osc(v9) < (5 (S^) and j osc{ip) {f,d^)<5{Yl) 

for some finite constant 6 (S^) < 00, i = 1,2, whose values do dot depend on the parameters 
(m,r/) G (/„ xV{En) and f G Osci 

Condition (-^2) is a rather basic and weak continuity type property. It states that the one 
step transformations of the flow of measures ()1.2p are weakly Lipschitz, in the sense that the 
mass variations and the integral differences w.r.t. some test function / can be controlled by 
the different initial masses and measures w.r.t. a collection of integrals of a possibly infinite 
number of test functions. It is satisfied for a large class of one step transformations r„. In 
section \2.'6\ we will verify that it is satisfied for the general class of Bernoulli and the PHD 
filters discussed in section II. 1.11 and section 11.1.21 

We are now in position to state the main results of this article. The first one is concerned 
with the exponential stability properties of the semigroup Tp,?! = 

(r^,,,r2„), withO <p<n 

associated with the one step transformations of the flow (|1.2p . A more precise description 
and the complete proof of the next theorem is provided in section [3l 

Theorem 1.1 We let ^p^n,^ ^''^d ^p^n,m semigroups associated with the one step 

transformations of the flow of total masses •— {•j^n-i) o,nd measures '■— 

r^(m„_i, .), with a fixed collection of measures v := {un)n>o S nn>o^(-^") '^^'^ masses 
m := (m„)„>o G Yln>o^n- When these semigroups are exponentially stable (in the sense 
that they forget exponentially fast their initial conditions) and when the pair of mappings 
Vn-i ^ ^n,u„-i '^^^ "^n-i ^ ^n.m„-i '^'^^ Sufficiently regular then we have the following 
contraction inequalities 

|ri^„(n',r?') -ri^„(n,r?)| V ||r2^„(n', r?') - T^^Ju, r/)||^^ < c e-^("-P) 

for any p < n, u,u' £ Ip, r],r]' G V{Ep), and some finite constants c < 00 and A > whose 
values do not depend on the time parameters p < n. 

The second theorem is concerned with estimating the approximation error associated 
with a A^-approximation model satisfying condition (jl.6p . The first part of the theorem is 
proved in section 11.21 The proof of the uniform estimates is discussed in section 13.11 (see for 
instance lemma [3^ . 



< c{n)\m-m'\+j \[r^ - l^l^^^, ^^,~^{d^) (1.14) 

< c{n)\m-m'\+[ \[v - v']{v)\ K,im',v')if, dip) (1.15) 
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Theorem 1.2 Under the assumptions (Hi) and {H2), the semigroup T^^^ satisfies the same 
Lipschitz type inequalities as those stated in Iil.l4\ ) o-n-d hl.l5\) for some collection of mea- 
sures and .) on B{Ep). In addition, for any N -approximation model satisfying 
condition lll.6\} we have the estimates: 

in 1 

E(|Fj'^(l)r)'<a.5^5(S^,J and E (|K''^(/)f ) " < a. 5 (E^ J (1.16) 

for any r > 1, and N > 1, with some constants Or < 00 whose values only depend on r. 
Furthermore, under the regularity conditions of theorem \1.1\ the couple of estimates stated 
above are uniform w.r.t. the time horizon; that is, we have that sup„>o Ylp=o ^ (^p,n) < 
for any i = 1,2. 

These rather abstract theorems apply to a general class of discrete generation measure- 
valued equations of the form (jl.ip . We illustrate the application of this pair of theorems 
in the analysis of the stability properties and the approximation convergence of the pair of 
multiple target filters presented in this introductory section. These results can basically be 
stated as follows: 

• The Bernoulli filter presented in section [1.1.11 with a sufficiently mixing prediction and 
almost equal survival and spontaneous births rates s„ ~ /^n(l) is exponentially stable. 

• The PHD filter presented in section 11.1.21 is exponentially stable for small clutter 
intensities and sufficiently high detection probability and spontaneous birth rates. 

• In both situations, the estimation error of any A^-approximation model satisfying con- 
dition ()1.6p does not accumulate over time. Furthermore, the uniform rates of conver- 
gence provided in theorem 11.21 allows to design stochastic algorithms with prescribed 
performance index at any time horizon. 

We end this section with some direct consequences of theorem 11.21 
Firstly, we observe that the mean error estimates stated in the above theorem clearly 
implies the almost sure convergence results 

lim Vnif) = Vnif) and lim 7^(/) = 7„(/) 

for any bounded function / on En- Furthermore, with some information on the constants 
Or, these L^-mean error bounds can be turned to exponential concentration inequalities. To 
be more precise, by lemma 7.3.3 in |5], the collection of constants in theorem [L2l can be 
chosen so that 

all < (2r-)! 2-"" /r\ and a^^+i < h^''+^{2r + 1)! 2-"" /r\ (1.17) 

for some b < 00, whose values do not depend on r. Using the above L.r-mean error bounds 
we can establish the following non asymptotic Gaussian tail estimates: 

F(|K-r?„](/)|>-^ + e) <exp(-^) with bn<bj2H^in) 
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The above result is a direct consequence of the fohowing observation 

Vr > 1 E (^J*^) < 6 ^ P (^7 > 6 + e) < exp {-e^/{2b)) 

for any non negative random variable U. To check this claim, we use the following Laplace 
estimate 

Vt > E (e**^) < exp (^^^ + bt^ ^F{U>b + e) < exp sup (^et - ^ 

It is worth noting that the above constructions allows us to consider with further work 
branching particle models in path spaces. These path space models arise in the analysis of 
the historical process associated with a branching models as well as the analysis of a filtering 
problem of the whole signal path given a series of observations. For instance, let us suppose 
that the Markov transitions M„ defined in (jl.lOl are the elementary transition of a Markov 
chain of the following form 

■■= (^p)o<p<„ ^En:= Yl K 

0<p<n 

In other words X„ represents the paths from the origin up to the current time of an auxiliary 
Markov chain taking values in some measurable state spaces E'^, with Markov transitions 
M^. We assume that the potential functions 5n,7„ only depend on the terminal state of the 
path, in the sense that (7„^^,^(X„) = g'^^^^^{X'^) , for some potential function g'n.'y,^ on E'^. In 
multiple target tracking problems, these path space models provide a way to estimate the 
conditional intensity of the path of a given target in a multi-target environment related to 
some likelihood function that only depends on the terminal state of the signal path. 

In practice, it is essential to observe that the mean field particle interpretations of these 
path space models simply consist of keeping track of the whole history of each particle. It 
can be shown that the resulting particle model can be interpreted as the genealogical tree 
model associated with a genetic type model (see for instance [5]). In this situation, rj^ is 
the occupation measure of a random genealogical tree, each particle represents the ancestral 
lines of the current individuals. 

We end this section with some standard notation used in the paper: 
We denote respectively by M.{E), V{E), and B{E), the set of all finite positive measures 
^ on some measurable space {E,£), the convex subset of all probability measures, and 
the Banach space of all bounded and measurable functions / equipped with the uniform 
norm ||/||. We denote by /~ and the infimum and the supremum of a function /. 
For measurable subsets A £ £, in various instances we slightly abuse notation and we 
denote f^{A) instead of /u(l^); and we set 6a the Dirac measure at a £ E. We recall that a 
bounded and positive integral operator Q from a measurable space {Ei,£i) into an auxiliary 
measurable space (£^2,^2) is an operator / 1— Q{f) from B{E2) into B{Ei) such that the 
functions 

x^Q{f){x) := [ Q{x,dy)f{y) 
JE2 

are <5i-measurable and bounded for some measures Q{x, .) G M.{E2). These operators also 
generate a dual operator ^ 1— )■ from M.[Ei) into M{E2) defined by {fiQ){f) := fJ-{Q{f))- 
A Markov kernel is a positive and bounded integral operator M with M(l) = 1. We denote 
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by Qp,n = Qp+iQp+2 ■ ■ ■ Qn, with p < n the semigroup associated with a given sequence of 
bounded and positive integral operator Qn from some measurable spaces {En-i,£n-i) into 
{En,£n)- For p = n, we use the convention Qn,n = Id, the identity operator. 

We associate with a bounded positive potential function G : x e E G{x) G [0,oo), 
the Bayes-Boltzmann-Gibbs transformations 

■■ r]&M{E)^^G{ri) eV{E) with ^ G(ri){dx) := G{x) r]{dx) 

provided r]{G) > 0. We recall that '^g{v) can be expressed in terms of a Markov transport 
equation 

vSr, = ^g{v) (1-18) 
for some selection type transition Sr){x,dy). For instance, we can take 

Sr,{x,dy) := ^ 5,{dy) + (^1 - ^) *(G-e)(r?)W (1.19) 

for any e > s.t. G{x) > e. Notice that for e = 0, we have S,^{x,dy) = "^Gi'n)idy). We can 
also choose 

Sr,{x,dy) := eG{x) S^{dy) + (1 - eG{x)) ^G{r]){dy) (1.20) 

for any e > that may depend on the current measure ry, and s.t. eG(x) < 1. For instance, 
we can choose 1/e to be the ?/-essential maximum of the potential function G. Finally, in 
the context of Bernoulli and PHD filtering we set JIn+i — A*n+i/Mn+i(l)) for any n > 0, the 
normalized spontaneous birth measures. 

2 Semigroup description 

2.1 The Bernoulli filter semigroup 

By construction, we notice that the mass process and the normalized measures are given by 
the rather simple recursive formulae 



^"^'^ ^ " (l-7n(l))+7n(l)^n(<7n) ^^"^^ ^ (1 - 7n(l)) + 7n(l)ryn(<7n) ^"^'^^ 

(2.1) 

and 

rin+i := an(7n) ^s„s„(??n)M„+i + (1 - a„(7„)) /I^+i 
with the mappings q;„ : 7 G M{En) 1-^ ck„(7) G [0, 1] defined by 

an(7) - 



l{Sn9n) + (1 - 7(l))/^n+l(l) 

By construction, if we set 7 = m x 77 then 
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with the collection of Markov transitions M„_|_i^^ defined below 

M„+i (7)M„+i(x,.) + (l-an (7)) (2.2) 

Next we provide an alternative interpretation of the mapping T^^^. Firstly, observe that 

^.„.„(^)M.+i,,(/) = ^§^±i4n| (2.3) 

rj iVn+l,m(,ljj 

with the integral operator 

Qn+l,m{f)ix) := mgn{x)Snix)Mn+lif){x) + (1 - m) Hn+lif) 

This implies that 

with the potential function 

Gn,m = mgnSn + - m) fin+l{l) (2.4) 

and the Markov transitions 

M ^9nSn /f^ , (l-m)^n+i(l) 

mgnSn + [l - m) IJ.n+l{l) mgnSn + - m) IJ.n+l{l) 

(2.5) 

The condition (Hi) is clearly not met for the Bernoulli filter (jl.Sp when s,i = and /i„+i(l) = 
0, since in this situation 7^ = for any n > 1. Nevertheless, this condition is met with 
In C (0, 1] and m0+_„(m) = 1, as long as s„ and ;U„+i(l) are uniformly bounded from below. 
It is also met for s„ = 0, as long as < /in+i(l) < 1 and the likelihood function given in 
(jl.Op is uniformly bounded. The condition is also met for Hn+i{l) = 0, as long as 7o(l) > 0, 
and the likelihood function given in (|1.9|) and the function s„ are uniformly lower bounded. 
We prove these assertions using the fact that 

7„+l(l) =7n(l) ^gArin){Sn) + //n+l(l) (2.6) 

with the updated mass parameters 7n(l) G [0, 1] given below 

^ 7n(l)??«(g«) 

(l-7n(l))+7n(l)??n(5n) 

If we set s~ := inf^;^ Sn and s+ = sup^^ s„ then 

Vn > 1 7„(1) G [m~,m+] 

with parameters 

"in = ^n(l) A s^_i and m+ = ^„(1) V (< 1) 

If Sn and /i„+i(l) are uniformly bounded from below then we have > 0. In addition, 
for the constant mapping s„, = ;U„+i(l), the total mass process is constant 

7„+i(l) = m++i = m-^-^ = /i„+i(l) 
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for any n > 0. Furthermore, in this situation the flow of normalized measures is given by 
the updating-prediction transformation defined by 

Vn>0 rin+i = ^^Mirjn)Mi^% 

(s) (s) 

with the UkeUhood function gn ' and the Markov transitions M^^i defined by 

9\t':=sn9n + {l-sn) and M^l/) -"^"^-^^(/) + (| " ^-j ^^-^(7) ^^.T) 

When Hn+iiX) = 0, the flow of normalized measures is again given by a simple updating- 
prediction equation 

rin+l = ^g^sSnn)Mn+l and 7n+l(l) = ^gn{'nn){Sn) X ^,7n(5n) (7n(l)) (2.8) 

with the increasing mappings Oa defined below 

X e [0, 1] ^ Oaix) := ax/[ax + (1 - x)] (2.9) 
In addition, if s~ > then 

> s„ X _ " > 

gnrUn + (1 - mn) 

as long as g~ := inf^;^ gn > 0, and 7o(l) > 0. We prove this inequality using the fact that 
the mapping (a,x) € [0,oo[x[0, 1] i— 6a{x) is increasing in both coordinates. In the case 
where s„ = 1, using the fact that and 6aO 6h = 6ab, we prove that 

7n+i(l) = e^^^g^) (7„(1)) = %^^o^^fe)(7o(l)) 

Conversely, when 7o(l) < 1 and < /n„+i(l) < 1 and s„ = 0, for any n > 0, then we 
have a constant flow of normalized measures 

Vn > 1 r]n='Pn 
and the total mass process is such that 

7n(l)e]0,l[ =^ 7n+l(l)=Mn+l(l)x [l - (7n(l))] e]0, 1[ 

with the convention JLq = rjo, for n = 0. In addition, if /li„_|_i(1) = 1 then we have 

72(n+i)(l) = %^=o(62p/&2p+i)(7o(l)) and 72n+l(l) = np=o'{&2p+i/62p)(7o(l)) 

for any n > 0, with the parameters bn ■= J^niOn)- We prove these formuae using the the fact 
that 1 — Oa{x) = 0i/„(l — x), and 9aO Of, = Oab- This again implies that m~ > as long as 
7o(l) > and the likelihood function are uniformly lower bounded. 
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2.2 The PHD filter semigroup 

By construction, if we set ^ = m x rj then we find that 

r^+i(?n,r?) = 7(5r„,^) +^„+i(l) and r^+i(m,r/) = ^-g,^ ^(r/)M„+i,^ 
In the above display, Mn+i^-y is the collection of Markov transitions defined below 



M„+i_^(3;, .) := On (7) •) + (1 - "n (7)) A^n.+i with a„(7) 



7(5n,7) 



7(5^,7) +/-in+l(l) 



The interpretation of the updating transformation ^ [rf) in terms of a Markov trans- 
port equation is non unique. For instance, using ()1.12p this Bolzmann-Gibbs transformation 
can be decomposed into two parts. The first one relates to the undetectable targets and the 
second is associated with non clutter observations. An alternative description is provided 
below. We consider a virtual auxiliary observation point c (corresponding to undetectable 
targets) and set 3^^ = 3^™ + ^c- We also denote by g^.^{.,y) the function defined below 

r„(l -dn) if y = c 

9l{-,y)={ ^ dngni.,y„) 



" K{y) + jidngni.,y)) 



if y / c 



In this notation, the updating transformation ^ {rj) can be rewritten in the following 
form 

= %,.,^(^) with = j y^{dy) gl{.,y) 

The averaged potential function 'g^^^ allows us to measure the likelihood of signal states w.r.t. 
the current observation measure y^. Using (jl.lSp . the Bolzmann-Gibbs transformation 
^ (^) can be interpreted as non linear Markov transport equation of the following form 

= ^-^"a ^"^{m^ri) = 7]Kn+i,^ with Kn+i,^ = Sn,-^Mn+i,.y (2.10) 

for some Markov transitions Sn,'y from En into itself. 

We also notice that condition {Hi) holds as long as the functions s„, 6„, and gni-,yn) are 
uniformly bounded and /i„,(l) > 0. It is also met when ^n(l) = 0, as long as r„ = (s„ + 6„) 
is uniformly lower bounded and 7^ or (i„ < 1. 



2.3 Lipsciiitz regularity properties 

Firstly, we mention that condition (-^2) can be replaced by the following regularity condition: 
{H2) : For any n > 1, f £ Osci(£'„), and any {m,r]), (m' ,r]') G (In x V{En)), the 
integral operators Qn,mri satisfy the following Lipschitz type inequalities: 

\\Qn,m7jif) - Qn,m'ri'{f)\\ <c{n) \m-'m'\+ j {[v - v']i'P)\ ^n,im' ,r,')if , dif) (2.11) 

for some collection of bounded measures ^n,(m' ,r]')U ■> •) ^{En) such that 

j OSc{(p) S„,(m,r,)(/, dtp) < 5 (S„) 

RR n° 7376 



Stability and Approximation of Branching Distribution Flows 



16 



for some finite constant 5 < oo, whose values do dot depend on the parameters (m, ij) G 
{InXV{En)) and f eOsci{En). 

We prove (-^2) ^ (|1.14p using the decompositions 

fT^^Qn,mr) fTl Tj Qn,m'ri' — TnTj ^Qn,mr} Qn,m'ri'~\ ~l~ ^rnTj 77T, 7^ j Qn,m'ri' 

and of course [mr] — m'r]'] = [m — m'] r} + m' [q — rj']. To prove (-^2) ^ (|l-15p . we let 7 = mr] 
and 7' = m'r]' and we use the decomposition 

[Vl{m, 7?) - Tl{m', V)] (/) = [iQnn - I'Qnn'] (/ " ^li^' . r}'){f)) 

The Bernoulli filter (jl.Sp satisfies {H'2)^ as long as the likelihood functions Qn given in 
(|1.9p are uniformly bounded above. In this situation, (|2.1ip is met with 

\\Qn,mr){f) - Qn,m'ri'{f)\\ < c{n) \m - m'\ + c' (n) \[r] - r]']{gn)\ 

for some finite constant c'{n) < 00. 

The PHD equation satisfies (^^2)5 long as the functions hn{y) + g'n.y with g'^ y := 
dngni-i y) are uniformly bounded above and below. To prove this claim, we simply use the 
fact that 



|5n,7 9n,j' II ^ Cji 



\m' — ml 



+ / ynidy) \[v'-riWn,y)\ 



This estimate is a direct consequence of the following one 

9n,^{x) -gn,yix) = / yn{dy) 



9n,y{x) - 7] {9n,y) 



Kiy) + l{9n,v) hn{y) + Y{9n 



n,yj 



Next, we provide a pivotal regularity property of the semigroup (rp,n)o<p<n associated 
with the one step transformations of the flow (jl.2p . 



Proposition 2.1 We assume that conditions (Hi) and {H2) are satisfied. Then, for any 
< p < n, f £ Osci(-E„), and any {m,r]), {m' ,r]') £ (Ip x V{Ep)), we have the following 
Lipschitz type inequalities: 



/ l['? 


-v']i'^)\ 


/ l[^ 


-v'm 



for some finite constants Cp{n) < 00, and some collection of bounded measures ^pn(m'r]') 
and ^l^n,{m',ri')(f' •) ^i^p) that 

osc((/7) < and J osc(^) ,^ (^_^)(/, < 5 (S^ „) (2.12) 

for some finite constant 6 (S^ < 00, i = 1,2, whose values do dot depend on the parame- 
ters {m,ri) G (Ip X V{Ep) and f G Osci 
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Proof: 

To prove this proposition, we use a backward induction on the parameter 1 < p < n. For 
p = {n — 1), we have r^_-^ ^ = F^, with i = 1,2, so that the desired result is satisfied for 
p = (n — 1). We further assume that the estimates hold at a given rank p < n. To prove 
the estimates at rank (p — 1), we recall that 

Fp_i,„(m, 7]) = Tp^n {Tp{m, rj)) ^ Vi = 1, 2 F^_^„(m, r/) = F^_„ (Fp(m, r])) 

Under the induction hypothesis 

|rp_i,„(rra,?7) -rp_i,„(m',r/')| = |F^ „ (Fp(m, r/)) - F^ „ (Fp(m', ry')) | 

< Cp{n) |F^(m, r/) - F^(m', r?')| 

+ 1 |[F2(m,r?)-F2(m',7?')]M| Sj,n,r,(„^'y)(«'^) 

On the other hand 
and 

I [rp("i, 1]) - Tl{m', 1]')] (if) I < c{p) \m-m'\ + j | [r? - r/'] | ((^, dV) 

The end of the proof is now clear. The analysis of ^p-i „ follows the same line of arguments 
and is omitted. This ends the proof of the proposition. ■ 



2.4 Proof of theorem 11.21 

This section is mainly concerned with the proof of the couple of estimates (|1.16p stated in 
theorem 11.21 

We use the decomposition 

(7^(l),r?;^)-(7„(l),r?n) = [Fo,„ (7o^(l), - Fo,„ (70(1), %)] 

n 

+ E [rp,n - rp-i,„ (7p^_i(l),ryp^_i)](2.13) 

p=i 

and the fact that 

Fp_i,p(7p^_i(l),7?p^_i) = (7p^(l),F2_,_p(7p^_i(l),r?p^_i)) 

to show that 

7^(i)-7n(i) = [rj,n (70^(1), %^)-rJ,„ (70(1),%)] 

p=l 
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Recalling that 7(^(1) = 7o(l)i using proposition 12.11 we find that 



n „ 

^/iV|7„^(l)-7n(l)| <E^p(") / ^Z'\d^) 



with the predictable measure 



p=0 



^1 



p,n,(in,ri) 



associated with the parameters (m, rj) 



(7p^(l),r2_, p (7p^_i(l),r?f_i)), with < p < n; and for p = 0, we set S, 



0,n 



^0,n,(7o(l),r?o)- 



Combing the generalized Minkowski's inequality with ()1.6p we have 



E 



for some constants Ur whose values only depend on the time parameter. This clearly implies 
that 

1 n 
p=0 

The normalized occupation measures can be analyzed in the same way using the decompo- 
sition given below: 



N 

Vn -Vn 



Kn (70^(1), %^)-rg,„(7o(l),%)] 

n 



p=l 



This ends the proof of the theorem 11.2 



3 Functional contraction inequalities 
3.1 Stability properties 

This section is concerned with the long time behavior of nonlinear measure- valued processes 
of the form (jl.2p . The complexity of these models depend in part on the interaction function 
between the flow of masses 7n(l) and the flow of probability measures r/„ = 7„/7„(l). One 
natural way to start the analysis of these models is to study the stability properties of the 
measure- valued semigroup associated with a fixed fiow of masses, and vice versa. These two 
mathematical objects are defined below. 

Definition 3.1 We associate with a flow of masses m = (m„)n>o £ nn>o-^" and probability 
measures v := {yn)n>o G Y\.n>o'^ i^ri) the pair of semigroups 

•= <I)i o o ^} and <I>^ •= <I>2 o o <t>? f3 1) 

p,n,u ■ n,Un-i ' ' ' l,uo "^i"-* ^p,n,m ' n,m„_i ' ' ' l,mo V"'-^/ 

with < p < n, and the one step transformations 

^n,v„_i : U e In-l ^ '^i,u„,iiu) ■■= {U, ^n-l) & In 
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By construction, using a simple induction on the time parameter n, we find that 
(mo,i^o) = (7o(l),??o) and Vn > 1 mn = ^i^^^_^{mn-i) and i^n = 

t 

Vn > (m„,z^„) = (7n(l),??n) 

In the cases that are of particular interest, the semigroups ^p^n,u and ^p^n,m ^i^l have a 
Feynman-Kac representation. These models are rather well understood. A brief review on 
their contraction properties is provided in section [321 Further details can be found in the 
monograph [S]. The first basic regularity property of these models which are needed is the 
following weak Lipschitz type property : 

{Lip{(^)) For any p < n, u,u' G Ip, r],r]' G ^{Ep) and f E Osci{En) the following 
Lipschitz inequalities 

< al,n (3-2) 

< <n/ \[V-V'm\ nln,^, if, dip) (3.3) 

for some finite constants Op „ < oo, with i = 1,2, and some collection of Markov transitions 
^pnrj' Z'^o"^ Osci{En) into Osci{Ep), with p < n, whose values only depend on the parame- 
ters p,n, resp. p,n and r]' . 

The semigroups ^p^n,u and ^p^n,m ™ay or may bot be asymptotically stable depending 
on whether Op „ tends to 0, as (n — p) — )• oo. In section 13.31 we provide a set of easily 
checked regularity conditions under which the semigroups associated with the Bernoulli 
models discussed in 12.11 are asymptotically stable. 

The second step in the study of the stability properties of the semigroups associated 
with the flow ()1.2p is the following continuity property: 

{Cont{^)) For any n > 1, u,u' £ In-i, ViV' ^ ^(-E'n-i) o-nd any f £ Osci{En) 

\^ir,iu)-'^i,'iu)\ < j \[n-n'\{V>)\^i,r,'{dip) (3.4) 

(/)| ^ \u-u'\ (3.5) 

for some finite constants < oo, with i = 1,2, and some collection probability measures 
^\iu' '^'^ Osci{En-i), whose values only depend on the parameters n, resp. n and v' . 

This elementary continuity condition allows us to enter the contraction properties of 
the semigroups ^],^n,v and ^'l^n,m the stability analysis of the flow of measures ()1.2p . 
The resulting functional contraction inequalities will be described in terms of the following 
collection of parameters. 

Definition 3.2 When the couple of conditions {Lip{^)) and {Cont(^)) stated above are 
satisfied, for any i = 1,2 and p < n we set 

'^p,n = ■^p+l ^p+l,n ^P,n = ^ O-p^q Og^n ^^"^ ^'p,n = ^ '^p,q '^\n (^-S) 

p<q<n p<q<n 



\%,n,M - %,n,M) 
\_^p,n,miji) '^p,n,m{^ )] (/) 
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The main result of this section is the following proposition. 

Proposition 3.3 // conditions {Lip{^)) and {Cont{^)) are satisfied, then for any p < n, 
u,u' G Ip, r],r]' G V{Ep) and f G Osci{En) we have the following Lipschitz inequalities 

|rp,„(n',??') -rp,n(^,??)| < cl^)i\u-u'\+cl'l j \[v-r]']if)\ ^l^n,u',ri'idf) 

|r^,nK,V)(/)-r^,n(n, '?)(/)] < c2;l |«-n'| 



1 









for some probability measures ^ ^, ^, (d^?) and Markov transitions ^pnm'r]'' 
lection of parameters 



c^'i = V c^'^ and c^'^ = V c^'^ 

'^p,n '^p,n ' / J P,g q,n u,ii,u, '-p,n / p,q q 



2,2 —1 

p<q<n p<q<n 



n—p 

^p,n — ^p,n + bp,ri Wk,rk+i 

1=1 p<ri<...ri<n l<k<l 
n—p 

'^pfn = ^p,n + Op,ri &rfc,rfe+i j with the convention r;+i = n. 

1=1 p<ri<...ri<n l<k<l 

In particular, the collection of parameters p < n introduced in /11.16\) and 

112.12]) are such that 

^ — '^p]n ^^'^ ^ (^p,n) ^ S!" 

The proof of this proposition is rather technical and it is postponed to section [531 in the 
appendix. Now we conclude this section with a direct application of the above estimates. 
The proof of the theorem 11.11 stated in the introduction and the uniform estimates discussed 
in theorem 11.21 are a direct consequence of the following lemma. 

Lemma 3.4 Suppose that r* = sup„>x < oo, and a^ ^ < Cj e~^^^'^~'P\ for any p < n, and 
some finite parameters Cj < oo and Aj > 0, with i = 1,2, satisfying the following condition 

Ai / A2 and ciC2 rV^ < (l - e'^^^^^^)) (e"(^^^^^) - e^^^^^^^)) 

Then, for any i,j G {1,2} we have 

( p(AiAA2) \ 

^ < c'' e"^(-^) with A = (Ai A X,) - log [l + crV^ ^_(,^^,^)_^_(,^^,^^ J > 
and the parameters d'^ defined below 



c 



- ci (l + c2'irV(e-^-e-^i)) c^'^ ^ ciC2T^ /{e~^ - e~^^) 



In particular, for any N- approximation models (7^ (1), r/^ ) of the flow (7„,(1), 77^) satisfying 
condition lil.b]) . the L^-mean error estimates presented in 111.16]) are uniform w.r.t. the time 
parameter 

\ 1 
supE(|yj'^(l)ry <ar ci'V(l-e~^) and supE f I V;^'^ (/)!") ^ < a, c^'Vll - e"^) 

n>0 ^ ^ n>0 ^ ^ 

wit/i some constants Or < 00 whose values only depend on r. 
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Proof: 

Under the premise of the lemma 

p<q<n p<q<n 

with c = C1C2 and r = r^r^. We further assume that Ai > A2 and we set A = |Ai — A2I. 

bp,n < cre-^^^^-^^-^^'+i)) Yl e-^(<'-(^'+i)) < cTe-^^«"-^)-(f+^)) /(I - e"^) 

p<(j<n 

In the same way, if A2 > Ai we have 

p<q<n 

This imphes that 

In much the same way, it can be shown that 

b'^^^ =< cr2e-(^i^^2)((n-i)-p)/(i _ g-A) (3 7) 

We are now in a position to estimate the parameters c^fn- Firstly, we observe that 

C^;^ < C2 e-^^("-f) + C2 ^ g-A,(n-p)g-(AiAA.)(n-n) 

1=1 \ / p<ri<...ri<n 

When Ai > A2, we find that 

«<..e-»'"-'g (^) ("7'') 

and therefore 

c^;^ < C2 e-^^(-^) (1 + cr Y^)"'" ^ ^p'-n = ^2 e-^("-^) 

with ^ 

A = A2-log (l + cr ^_J_'^_,^ ) >0 

as long as 

cr < (^1 - e-^^) _ g-Ai~J 

When A2 > Ai we have A2 = Ai + A, we find that 

"~P / 2Ai \ ' 

c^;^<C2e-^^("--)+ 026-^^(^-^)5: (^^) E -"^^"'"^ 

(=1 p<ri<...ri<n 
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from which it fohows that 



Using a similar hne of argument as above, we have 



p,n 

with 



A = Ai - log I 1 + cr — ^ I > 

as long as 



g-Ai _ g-A2 



We conclude that 
with 



/ e(AiAA2) \ 

A = (Ai A A2) - log (^l + cr ^_(^^^^^^_^_(^^^^^^ j >0 

as long as 

CT- < ^1 - e-(^i^^2)^ ^g-(AiAA2) _ g-(AiVA2)j 

Using (j3.7p we also show that 

c^'i < c^'i e-^("-^') with c2'i = cr2 



P'" - e-(AiAA2) _ e-(AiVA2) 

Using these estimates 

p<q<n 

and 

p<g<n 

Since Ai > A we find that 

< ci e-^i("-P) + c^'iciT^ e-^""-!)-^')/!! " e"^') with A' = Ai - A > 
This yields 

4;^ < c^'^ e-"("-P) with ci'^:=ci (l + c^' V/(e-" - e-^^; 
Finally, we observe that 

p<q<n 

which implies that 

cl'^n < c^'^ e-^("-P) with := cTV(e-^ - e'^i) 
This ends the proof of the lemma. 
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3.2 Feynman-Kac models 

We let Qp^m with < p < n, be the Feynman-Kac semi-group associated with a sequence of 
bounded and positive integral operator Q„ from some measurable spaces {En-i,£n-i) into 
{En,£n)- For any n > 1, we denote by and M„ the potential function on En-i and 

the Markov transition from En-i into En defined below 

Gn-i{x) = Qn{l){x) and M„(/)(x) - ^"^^^^^^ 



Qn{l){x) 



We also denote by ^p^m < p < n, the nonlinear semigroup from V{Ep) into V{En) defined 
below 

Vt? G V[Ep), V/ G B[En) %,nmf) = VQp,n{f)/vQp,na) (3-8) 

As usual we use the convention ^n,n = Id, for p = n. It is important to observe that this 
semigroup is alternatively defined by the formulae 

^P,n(r/)(/) = with Gp,„ = gp,„(l) and Pp,„(/n) = Qp,„(/n)/Qp,n(l) 

The next two parameters 

rp,n= sup ^^''^j^}. and /3(Pp^„) = sup \\Pp,n{xp, ■) - Pp,n{yp, ■)\\tv (3.9) 

measure respectively the relative oscillations of the potential functions Gp^n and the contrac- 
tion properties of the Markov transition Pp^n- Various estimates in the forthcoming sections 
will be expressed in terms of these parameters. For instance and for further use in several 
places in this article, we have the following Lipschitz regularity property. 

Proposition 3.5 (|6|) For any fn G Osci(£'„) we have 

li^pAVp) - %,n{l^p)] ifn)\ < 2 rp,n P{Pp,n) | " /"p] ^'^^^(/n) | (3.10) 

for some function Pp\{fn) G Osci{Ep) that doesn't depends on the measure r]p. 

Our next objective is to estimate the the contraction coefficients rp^n and l3{PpA in terms 
of the mixing type properties of the semigroup 

Mp^Xp, dXn) ■■= Mp+iMp+2 ■ ■ ■ Mn{Xp, dXn) 

associated with the Markov operators M„. We introduce the following regularity condition. 



{MG)m There exists an integer m > 1 and a sequence (ep(Af ))p>o G (0, 1) and some 
finite constant rp such that for any p>0 and any {x,x') G Ep we have 

Mp^p+mixp, .) > ep{m) Mp^p+rnix'p, .) and Gp{x) < rp Gn{x') (3.11) 



It is well known that the above condition is satisfied for any aperiodic and irreducible 
Markov chains on finite spaces. Loosely speaking, for non compact spaces this condition is 
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related to the tails of the transition distributions on the boundaries of the state space. For 
instance, let us suppose that En = M and M„ is the bi-Laplace transition given by 

Mn{x,dy) = ^ e-c{n)|j/-A„(x)l 

for some c(n) > and some drift function An with bounded oscillations osc(^„) < oo. In 
this case, it is readily checked that condition (M)m holds true for m = 1 with the parameter 
en-i(l) = exp(-c(n) osc(^n))- 

Under the mixing type condition {M)m we have for any n > m > 1, and p > 1 

rp,p+n < ep{my^ Vp+k (3.12) 
0<A:<m 

and 

[ra / m] — 1 

f^{Pp,P+n)< n (l-4+U ^ith 6f):=6^M n (3.13) 

k=0 0<k<m 

Notice that these estimates are also valid for any n > 0. Several contraction inequalities 
can be deduced from these estimates (see for instance chapter 4 of the book [5]). To give a 
flavor of these results, we further assume that {M)m is satisfied with m = 1, and we have 
e = inf„ e„(l) > 0. In this case, we can show that 

I p,p+n <rp/e and /3(Pp,p+„) < (l - e^)" 
We end this short section with a direct consequence of proposition 13. 5[ 

Corollary 3.6 Consider the Bernoulli semigroup presented in section \2.1\ For constant 
mappings Sn = /in+i(l); the first component mapping is constant ^^_^_^j^^{u) = Sn and the 



seconc 



d component mapping ^n+i m„iv) — '^gM(.v)^n+i induces a Feynman-Kac semigroup 

with the likelihood function gn '' and the Markov transitions M^^^ defined in \2. 7\ ). In this 
situation, the condition iS.S^) is clearly met with a^ ^ = 0, for anyp < n,. We further assume 
that the semigroup of associated with the Markov transitions Mn satisfies the mixing property 
stated in the l.h.s. of 13.11\) for some integer m > 1 and some parameter ep{m) g]0, 1]. In 
this situation, the condition i3. 3|) is also met with the collection of parameters a^ n given 
below 

[(n— p)/mj — 1 

ain < 2 pm) n (i - 

with 

Pp{m) := e^^ {m) J| r^(sfe)rfc(l) and e^p''""''^ = el{m)rp{sp) / rk{skfrk{l) 

p<k<p+m p<k<p+m 

and the collection of parameters rn{sn) defined below 

I \ ^n9n ~l" (1 ~ ^n) /^\\ 

rn{Sn) ■■= . T (< rn[l)) 

SnQn + V-L Sn) 
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3.3 Bernoulli models 

This section is concerned with the contraction properties of the semigroups ^p^n,^ ^p^n,m 
associated with the Bernoulh filter discussed in section [2m Before proceeding, we provide 
a brief discussion on the osciUations of the hkehhood functions gn given below 

gniXn) = (1 - dniXn)) + dn{Xn)yn {ln{Xn, O/^n) 

in terms of some [0, l]-valued detection probability functions d„, some local likelihood func- 
tions /„, and some positive clutter intensity function The oscillations of these likelihood 
functions strongly depend on the nature of the functions In). 
Assuming that > we have 

(1 - + dr ^ yn (1) <g-<g^<{l- + C'+ ^ yn (1) (3.14) 

with the parameters 

dn"^ = dt ^i+yr,{i)>h- + dn l/+y„(i)<ft- 

dn = dn li-y„(l)>/i+ + dt l/-y„(l)</i+ 

The semigroup contraction inequalities developed in this section will be expressed in terms 
of the following parameters 

5„(5<7):=44, '^n(<?):=4 and 8'n{g) := ^ ^ gt^ 

Qn Sn Qn 9n 

For time homogeneous models {dn,hn,ln) = id,h,l), with constant detection probability 
dnix) = d and uniformly bounded number of observations sup^3^„(l) < 3^^(1) < oo we 
have the following estimates 

il-d)< g- <g+ <{l-d) + d^y+ (1) 
In this situation, we have 

^ ' + 1^ f 

For small clutter intensity function with > and l^ > we also have the observation 
free estimates % < , from which we find that the upper bound 



n>0 

and for d < 1 



6(9) := sup 6nig) ^ 1 + ^ 3^ (1) ' (3-15) 



6'ig) :=sup5;(<7) <sup|(l-(i) + d-^ 3^(1),^| (3.16) 



To be more precise, if we set inf„3^„(l) = y (1) then 

1 < ^yi^y ^ ^'(9) <ii-d)+d^y+ (1) 
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In addition, if we have d{l — d)y{l) < h and ti < 1 then we find the the observation 
free estimates 

dy{l) l+/h- < 1/(1 - d) ^ 5'{g) <{l-d) + ^ 
Conversely, we have the observation free estimates 

^^y(l)+ < 1 ^ 5'{g) < -, < 

We are now in position to state the main result of this section. 

Theorem 3.7 // //„+i(l) g]0, 1[, < s~ < s+ < 1, and the semigroup Mp^n satisfies the 
condition stated in the l.h.s. of h3.11\) for some integer m > 1 and some positive constant 
ep{m), then the condition {Lip{^)) is met with 

[n/mj — 1 

ap,„ < 2 ep^6'p{g) (l - el) and „ < 2 pp{m) JJ (^1 - e^^l^ 

p<k<p+n k=0 

with some parameters 

s- /i„+i(l) 1-4 l-^„+i(l) 



En > inf 
and 



/i„+l(l)' st 'l-Mn+l(l)' 1 



Pp{m) < ep{m) ^ 5p+k{sgf and e^"") > ep{rrif 6p{sg) ^ Sp+kisg)' 

0<k<m 0<k<m 

In addition condition (Coni($)) is met with 

Tn+l<6n{g) [{S^ - S~) + \\Sn - tJ-n+l{l)\\] and T^+i < (5^ (5) SUp | 



Sn '/^n+l(l) 

The proof of the theorem is postponed to section 15.21 To give a flavour of these estimates 
we examine time homogeneous models 

with constant detection and survival probabilities dn{x) = d, Sn(x) = s, and uniformly 
bounded number of observations sup„3^n(l) < 3^(1) < oo. In this situation, we have 
{ep{m) , ep^\m)) = {e{m) , e^^^m)) and using the estimates p.l5p we prove the following 
bounds 

r^+i < S{g) \s - /.(l)l and r^^^ < 6'{g) 

and 

al^<2e-'5'{g){l-ey and < 2e{m)-'5{gf- (l - 6(^)^5(5)-^-+^) ^"/"^J 
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with some parameter e such that 

^ Ml) 1-^ ^ 

Ml)' s 'i-Mi)' /-'- 

It is also readily verified that the assumptions of lemma [33] are satisfied with the parameters 

< 5{g)\s-^{l)\ r2 < 6' (g) {{(,{1) V s) / A s)) 

ci = 2e-i(5'(c/) C2 = 2e(m)-i (1 - e(m)25(c/)~5™+i) ^ S^g)'^"" 

and the Lyapunov constants 

Ai = -log(l -e^) and A2 = -- log (l - e(m)2%)-5'"+i) 

m 

We notice that e tends to 1 and tends to 0, as |s — /u(l)| tends to 0. Thus, there exists 
some <f > such that 

Ai > A2 and ciC2T^t'^ < ~ ^~^') (^~^' ~ ^"^0 

as long as |,s — /i(l)| < We summarize this discussion with the following corollary. 

Corollary 3.8 Consider the time homogeneous model discussed above. Under the assump- 
tions of theorem \3. 7\ for any N -approximation models (7,^(1),??^) of the Bernoulli model 
(7„(l),77n) satisfying condition lll.6\) . the hr-mean error estimates presented in hi. 10(1 are 
uniform w.r.t. the time parameter 

supE('|yj'^(l)r)^ <ar ci'V(l-e-^) and supE f |v;^'^(/)r) ^ < a, c^'VCl - e"^) 

n>0 ^ ' n>0 ^ ' 

with the parameters (c^'^, c^'^, A) defined in lemma and some finite constants Or < 00 
whose values only depend on r. 

Remark 3.9 When ;U„+i(l) = u>e have seen in 112. 8\) that 

^n+i,^.„H = ^9nK)(sn) X and $2^i_^Jr/) = ^'g„,„(r/)M„+i 
with the collection of mappings 9a, with a G [0, oo[, defined in 112.9]) . Using the fact that 
K+i u in) - $1+1 , = "^aM^Sn) V {g ) _ ^ 

1 n+l,u^y ) n+M )\ + (1 _ ^)] + (1 _ ^/)] I I 

one proves that \3. is met with the rather crude upper bound 

ap,„ < Yl "'k,k+i and 4,fc+i < (4^^ )/(l A 5fe )^ 

p<k<n 

We also notice that the second component mapping ^^+1 m„ doesn't depends on the param- 
eter nin, and it induces a Feynman-Kac semigroup of the same form as the one studied in 
section Assuming that the mixing condition stated in the l.h.s. of is satisfied 

some integer m > 1 and some parameter €p{m) > 0, one can prove that 113.!^) is met with 
the collection of parameters a^ ^ given below 

[(n— p)/mj — 1 

4,n < 2 Pp(m) Yl (l-4+L) ^^^^ Pp{m) = ep\m) JJ 5g{sg) 

k=0 p<q<p+m 

and the collection of parameters e^"*^ = e^"*^ = ^'^p{m) / Wp^^q^pj^m^Q^^d) • 
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3.4 PHD Models 

This section is concerned with the contraction properties of the semigroups ^p^n,^ ^p,n,m 
associated with the PHD filter discussed in section 11.1.21 and in section 12.21 

The analysis of these nonlinear models is much more involved than the one of the 
Bernoulli models. We simplify the analysis and we further assume that the clutter intensity 
function, the detectability rate as well as the survival and the spawning rates introduced in 
section 12.21 are time homogeneous and constants functions, and we set 

(6„(x),/i„(x),s„(x),r„(x)) = (6, /i,s,r) 

To simplify the presentation, we also assume that the state spaces, the Markov transitions 
of the targets, the likelihood functions and the spontaneous birth measures are time ho- 
mogeneous, that is we have that En = E, E^ = E^ , = M, gn{x,y) = g{x,y) and 
Hn+i = fJ-- Without further mention, we suppose that r(l — d) < 1, /x(l) > 0, r > 0, and for 
any y S E^ we have 

< g~{y) := ini g{x,y) < g+{y) := sup g{x,y) < oo 

Given a mapping 6 from E^ into M, we set y~{6) ■= inf„3^„(0) and y~^{9) := sup„3^n(^)- 
We recall from (jl.lOp that the PDH filter is defined by the measure-valued equation 

7n+l = 7nQn+l,j„ 

with the integral operator 

l) = anau{Xn)Mn+l{Xn-,dXn+l) +7n(l) ^^n+l{dXn+l) 

with the function gn,'y„ defined below 

9{x,y) 



9n,f„{x) = r{l - d) + rd / 3^„(rfy) 



h + d'yn{g{.,y)) 



We also notice that the total mass process and the normalized distribution fiow are given 
by the following equations 

7n+l(l) = ^n+l,^„(7n(l)) 

= 7n(l) r{l-d) + j yn{dy) u;^„(i)(r/„, y) + ^(1) 

DC 7n(l) r{l-d) r]nM + j 3^„((iy) w^„(i)(?7„, y) ^'g(.^j^)(?7„)M + /i(l) JI 
with the probability measure Ji and weight functions w defined below 
]l{dx) = nidx) / IJ.{1) and y) := r ( 1 



h + duri{g{.,y)) 
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For null clutter parameter h = 0, we already observe that the total mass transformation 
^n+i,ri„ doesn't depend on the flow of probability measures r/„ and it is simply given by 

^n+i,,„(7n(l)) = 7„(1) r(l -d)+r 3^„(1) + /i(l) 
In this particular situation, we have 

7^:^(1) =7n(l) = (r(l-d)r7o(l)+ E {r{l-d)r-'-Hryk{l)+Kl)) 

0<k<n 

Now, we easily show that the pair of conditions ()3.2p and (j3.4p are satisfied with the param- 
eters Op „ = (r(l — d))'^~^ and = 0. In more general situations, the total mass process is 
not explicitly known. Some useful estimates are provided by the following lemma. 

Lemma 3.10 We assume that the number of observations is uniformly bounded; that is, we 
have that 3^^(1) < oo. In this situation, the total mass process 7n(l) and any approximation 
model 7^(1) given by the recursion (with the initial condition 7o^(l) = 7o(l)J take 

values in a sequence of compact sets In C [m~ ,m~^] with 

~i n j\ [ ^ + rd y and m/ := 7o(l) + 



l-r(l-d) V \h + d^i{l)g- J J ' '^^ ' I - r{l - d) 

Proof: 

Using the fact that 7ra(l) > ^(1) we prove that 

"-['-h + dJ^D ,-(,)) 

from which we conclude that 

7„(1) r(l -d)+r yh,n{l) + < <I>i+i,,„(7n(l)) < ln{l) r{l - d) + r 3^„(1) + ^(1) 
with the random measures 

dfj'i^) g~{y) 



yhAdy) ■■= yn{dy) 



h + dfi{l) g-{y) 



For any sequence of probability measures u := (fn)n>o G V{E)^, and any starting mass 
u £ [0, c«[ one conclude that 

f n J. ry^{l)+K'^) ^ ^1 r n ^ ry+{l) + fi{l) 

(r(l-d)) u+ < ^o,n..N < - rf)) n+ 

This implies that 7^(1), 7^(1) € In C [m~,m+] with 

m := — — = — \ l + rdy ' 



l-r(l-d) l-r(l-d) V \h + dfi{l)g- 

The end of the proof of the lemma is now completed. 
We are now in position to state the main result of this section. 
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Theorem 3.11 We assume that the number of observations is uniformly bounded; that is, 
we have that < oo. In this situation, the condition {Lip{^)) is met with the Lipschitz 

constants a^^^^ < Ylp<k<n'^k k+i' ^^^'^ ^ = and the sequence of parameters ("ran+i)„>o' 
i = 1,2, defined below 



and 



a^,„+i < r(l -d)+ rdh J„ 

m) [{^-d) + dyJj,^£^ f 



[h + dm g 

+ hdyn 



9^-9 
{h+din~ g~)'^ 



(1 -d)m-+ dm-y^ (Tl+i^ ) + ^(1)/ 
In addition, condition {Cont{^)) is met with the sequence of parameters 

[i-d) + hd yn ' — 



< rdhm'^ yn 



9 -9 



[h + dm g 



{h+dm~ g~)'^ 



(.. 



{l-d)m +dm yn y h+dm-g- 
The proof of theorem 13.111 is postponed to section 15.41 



+ /i(l)/r 



Corollary 3.12 We assume that y^ {9^/9^) md y^ (5^/(5 )^) < oo- In- this situation, 
there exists some parameters {) < kq < \, ni < oo, and K2 > Q such that for any d > kq, 
fi{l) > Ki, and h < K2, the semigroups ^p^n,u ^''^d ^p^n,m satisfy the pair of conditions 
{Lip{^)) and {Cont{^)) with some parameters (Op^n) '^n)j=i,2,p<n, satisfying the assumptions 
of lemma In particular, for any N -approximation models {'yni^)i'nn) of the PHD 

equation (7n(l))??n) satisfying condition U.6\) . the hr-mean error estimates presented in 
U.16\) are uniform w.r.t. the time parameter 



supE(|yj'^(l)|'y <ar ci'V(l 

n.>0 



and supE ( ) < a, c^'Vll - e"'') 

n>0 



1,2 ^2,2 



A) defined in lemma and some finite constants Ur < 00 



with the parameters (c 
whose values only depend on r. 

Proof: 

There is no loss of generality to assume that r(l — d) < 1/2 < d and > I > h. Recalling 
that m~ > /u(l), one readily proves that 



m ' 



7o(l) 



+ 



1 



1 + 



VI , 



Ml) Ml) l-r(l-d) 
If we set d{g) := pVy+ (^|^) V y+ 

2h 

ai^n+i/r <{l-d) + 6{g) and a^,„+i/r < 

as well as 

2h 1 

Tn+i/r < S{gf and r^_^i/r < — 



3^+(l) <2 + 7o(l) + 2rJ^+(l) :=p 



then we find the rather crude estimates 

2h + l3{M) 



/3(M)(1 -d) + 



{l-d) + 



2h 



Si9) 
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from which we find that 



rV2< 



2hr^ 



{l-d) + 



2h 



Ha) 



Ho? 



(3.17) 



Thus, there exists some ^ < hq < 1 and some ki < oo so that for any d > kq and any 
/i(l) > Ki we have 



a, 



1 

n,?i+l 



< r 



< r 



(1 - d) + 
(1 - d) + 



Mi)2 

3 



5(g) := e-^i < 1 

J(5) := < 1 ^i^j^ < A2 < Ai 



Finally, using (|3.17|) we find some K2 > such that for any /i < K2, we have that r^r^ < 
(1 - e-^^) (e-^2 _ g-^i). The end of the proof is now a direct consequence of lemma [37 
This ends the proof of the corollary. 



4 Stochastic particle approximations 

4.1 Mean field interacting particle systems 
4.1.1 Description of the models 

The mean field type interacting particle system associated with the equation (|1.2p relies on 
the fact that the one step mappings T^^^^ can be rewritten in the following form 

r^+l(7n(l),f?n) = ??ni^n+l,7„ with 7„ = 7n(l) X (4.1) 

for some collection of Markov kernels Kn+i^-^ indexed by the time parameter n and the set 
of measures 7 G Mj^{En). We mention that the choice of the Markov transitions Kn^'y is not 
unique. In the literature on mean field particle models, Kn.-y are called a choice of McKean 
transitions. Some McKean interpretation models of the Bernoulli and the PHD filter models 
p. 81 ) and (|1.1U|) are discussed in section [22] (see for instance (j2.10p ) and in section [XT] (see 
for instance 12. 2p 

These models provide a natural interpretation of the distribution laws r]n as the laws 
of a, non liiiGBi ]VIa,rkov cha;in whose clenienta^ry transitions ^ Xn+i depends on 
the distribution = Law(X„), as well as on the current mass process 7n(l)- In contrast 
to traditional McKean model, the dependency on the mass process induce a dependency 
of all the fiow of measures 77^, for < p < n. For a thorough description of these discrete 
generation and non linear McKean type models, we refer the reader to [S]. 

In further developments of the article, we always assume that the mappings 

are pointwise known, and of course measurable w.r.t. the corresponding product sigma 
fields, for any n > 0, > 1, A^+i £ £n+i, and any Xn E En- In this situation, the mean 
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field particle interpretation of this nonlinear measure- valued model is an -valued Markov 
chain f!^^ = , with elementary transitions defined as 

V /l<i<N 



7;Vi(l) = 7/,'(l)<(G„,,„.) (4.2) 

N 

p(cl^lGdx 14"^)) = n ^-+i,7-(ef'^^^*) (4-3) 



i=l 



with the pair of occupation measures {jn ^Vn) defined below 

1 ^ 

i=i 

In the above displayed formula, stands for the cr-field generated by the random sequence 
{ip^^)o<p<n, and dx = dx^ x . . . x dx^ stands for an infinitesimal neighborhood of a point 
X = {x^ , . . . ,x^) G . The initial system ^^"^ consists of N independent and identically 
distributed random variables with common law r/Q. As usual, to simplify the presentation, 
when there is no possible confusion we suppress the parameter N , so that we write and 
instead of ^1^^ and ^n^'*'' . 

4.1.2 Convergence analysis 

The rationale behind the mean field particle model described in (|4.3p is that ^/^^ is the em- 
pirical measure associated with N independent variables with distributions K^j^-^^^n (^J^, dx) , 
so as long as 7^ is a good approximation of 7„ then ?7^x should be a good approximation 
of Roughly speaking, this induction argument shows that rj^ tends to r/n, as the 

population size N tends to infinity. 

These stochastic particle algorithms can be thought of in various ways: From the physical 
view point, they can be seen as microscopic particle interpretations of physical nonlinear 
measure-valued equations. From the pure mathematical point of view, they can also be 
interpreted as natural stochastic linearizations of nonlinear evolution semigroups. From the 
probabilistic point of view, they can be interpreted as a interacting recycling acceptance- 
rejection sampling techniques. In this case, they can be seen as a sequential and interacting 
importance sampling technique. 

By construction, the local fluctuation random fields {W^)n>Q defined in (jl.Sp can be 
rewritten as follows 

Using Khintchine's inequality, we can check that (jl.6p is met for any r > 1 and any /„ G 
Osci with the collection of universal constants given below 



< (2r)! 2-'' /r\ and a^^+i < (2r + 1)! 2''' /r\ 

We end this section with a brief discussion on the PHD equation presented in (jl.lOp . 
This model combines in a single step the traditional updating and a prediction filtering 
transition. This combination allows us to reduce the fluctuations of the local sampling 
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errors and their propagations w.r.t. the time parameter. Since these updating-prediction 
models are often used in the hterature of multiple target tracking, we provide below a short 
summary. If we set 



' (1 - d„) if y = c 

if y^c 



, hn{y)+i{dngn{-,y)) 

then 

7„+i = 7„Qn+i + with Q„+i(/) := r„ M„+i(/) 

with the updated measures defined below 

7n(/) := 7n(9n,7„/) ^ith = j y^{dy) gn,^„{.,y) 

Notice that 

7n(l) = 7n(5n,7„/) and 9jn{dx) := %{dx)/%{l) = '^gc^ Jr]n){dx) 
from which we find the recursive formulae 

7n(l) A updating / ^^(i) \ prediction / 



Vn J \ Vn J V ^n+l 

with the prediction transition described below 

7n+i(l) = %{rn) + ^n+i(l) and rjn+i = iVn) M'„+i,% 
In the above displayed formula, M^_^-^ is the Markov transition defined by 

^n+l,%i^' •) = «n(7n) Mn+l{x, •) + (l " a'niln)) 

with the collection of [0, l]-valued parameters 0^(7^) = 7n(?'n)/(7n(?'n) + /^n+i(l))- It 
should be clear that the updating and the prediction transitions can be approximated using 
a genetic type selection and mutation transition. Each of these sampling transitions intro- 
duces a separate local sampling fluctuation error. The stochastic analysis of the correspond- 
ing mean field particle interpretations can be developed using the same line of arguments 
as those used for the particle model discussed above. 

4.2 Interacting particle association systems 
4.2.1 Description of the models 

We let {An)n>o be a sequence of finite sets equipped with some finite positive measures 
{t^n)n>o- We further assume that the initial distribution 70 and the integral operators 
(5n+i,7„ in (|l-ip have the following form 



70 



{da) qSi,^^ 
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In the above display tjq stands for a collection of measures on Eq, indexed by the parameter 

a G and Q^+i -y^ is a collection of integral operators indexed by the parameter o E ^n+i- 
In this situation, we observe that 

7o(l) = fo(l) and rjQ = j AQ{da) rif' with Ao{da) := vo{da) / 

We also assume that the following property is met 

Qiti,7(l) « Ql"|i,,(/)/gSi,,(l):=M("^\(/) (4.4) 

for some function Gn"^ on and some Markov transitions M^^^ from into E^+i whose 
values do not depend on the measures 7. For clarity of presentation, sometimes we write 
^'1^^ instead of ^'^(a) . 

tin Gji 

Definition 4.1 We consider the collection of probability measures r]li"^ G V{En), indexed 
by sequences of parameters 

a„ = (ao, . . . , an) G -4[o,„] := (A x . . . x An) 

and defined by the following equations 

ryi-) = ($(-)o...ocI>(-))(r^('^o)) (4.5) 

with the mappings : V{En-i) — )• V{En) indexed by a £ An and defined by the updating- 
prediction transformation 

We illustrate these abstract conditions in the context of the multiple target tracking equation 
presented in p.iop . In this situation, it is convenient to add a pair of virtual observation 
states c, c' to E^ . Using this notation, the above conditions are satisfied with the finite sets 
An+i and their counting measures Un+i defined below 

An+l = {Yll<i< Nl} U {C, C'} Vn+l =yn + 5, + G M{An+l) 

Using (jl.lOp and ()1.12p . we check that ()4.4p is met with the couple of potential functions 
and Markov transitions defined by 

{{rndngn{-,yn),Mn+l) for 2/„ {c, c'} 
(r„(l-d„),M„+i) for yn = c 
(l,7^„+i) for Vn = c' 

In this case, we observe that 

Q'^n+l^'y^i^n.) •) = G^^^^{Xn) M^_^^{Xn, •) 
(v ) 

with the potential function Gri,7„ defined below 

{[hniVn) + ln{dngn{-,yn))]~^ for y„ {c, c'} 
1 for yn = c (4.6) 

/u„+i(l)/7„(l) for y„ = c' 

Under our assumptions, using (jl.2p . we have the following result. 
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Proposition 4.2 The solution the equation il.^} has the following form 

7]n = / An{da) T]'^^ 



with a total mass process 7n(l) and the association measures G 'Pi-^[o,n]) defined by the 
following recursive equations 

7„+i(l) = 7„(1) r]niGn,^„) and An+i = ^n+i (7„(1), A„) 
With the mapping 

Un+i : {m,A) £ (]0,oo[xP(^[o,n])) ^ l^„+i(m, A) G PM[o,n+i]) 
defined by the following formula 

a„+i im,A) {dia,b)) Aida) r/i") (g^'^ ,^ (4.7) 

Proof: 

The proof of the above assertion is simply based on the fact that 



r]n+i oc / = / Anida) Vn+i{db) 



This clearly implies that 

Tl [m,j A{da) r^^:\)) = j ^n{m, A) {d{a,b)) vt'^ 
This ends the proof of the proposition. ■ 

By construction, we notice that for any discrete measure A G 'Pi-^[o,n-i])i and any 
collection of measures r/(") G V{En-i), with a G ^[o^„_i] we have the formula 

Tl [m,j A{da) = j {m, A) id{a,b)) cl>(f) (ry('^) 

4.2.2 Particle approximation models 

To get some feasible solution, we further assume that rjn^ ^G^n^7n^ are explicitly known for 

any sequence of parameters (o, b) G (-^[o,n] ^ -^n+i) ■ This rather strong condition is satisfied 
for the multiple target tracking model discussed above as long as the quantities 

vk^°''""-'Hrndngn{.,yn)) r?(«0,?/0,...,?/n-i) (^^(i _ rf^)) r?(»« '^0'-'^"-) (d^^n ( • , ^n)) 

are explicitly known. This condition is clearly met for linear gaussian target evolution and 
observation sensors as long as the survival and detection probabilities Sn and dn are state in- 
dependent, and spontaneous birth /I^ and spawned targets branching rates bn are Gaussian 
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mixtures. In this situation, the collection of measures ^^("O'^O'---'^" are gaussian distribu- 



tions and the equation (j4.5p coincides with the traditional updating-prediction transitions 
of the discrete generation Kalman-Bucy filter. 

We let Aq = jj "^^^i Sgi^ , be the empirical measure associated with N independent 
and identically distributed random variables (ag)i<j<Ar with common distribution Aq. By 
construction, we have 



with some local sampling random fields satisfying (jl.6p . We further assume that 7o(l) is 
known and we set 7,^ = 70(1) tjq . 

7f (1) = 7o"^(l) V^iG,^,^) and ryf := J Af {da) r^^^^ 

with the occupation measure A^ = 'Yld=i associated with N conditionally inde- 
pendent and identically distributed random variables a,\ := (aQ^,a^^) with common law 
r^i (7q^ (1),^q^). By construction, we also have 

■= I ^1 (7o"^(l),0 (da) + < = T? {l^ilU^) + < 



with some local sampling random fields satisfying ()1.6p . Iterating this procedure, we define 
by induction a sequence of A^-particle approximation measures 

7n^(l) = 7^-1(1) , ) and r?^ := / A^ (da) r^^^ 



with the occupation measure A^ = Yli=i ^a-i^ associated with N conditionally independent 
and identically distributed random variables a^ := (og a^^ „) with common law 

(7^1(1), • Arguing as above, we find that 

r/^ = /^^n (7:-i(l), All) {da) r?^ + ^ = Vl (7„^_i(l), + I^.^ 



with some local sampling random fields satisfying ()1.6p . 
4.2.3 Convergence analysis 

The main objective of this section is to show that A^-particle occupation measures A^ 
converge in a sense to be given, as tends to 00, to the association probability measures 
An- To this end we observe that the one step mapping ^n+i introduced in (|4.7p can be 
rewritten in the following form 

O A\(T?\ ^Qn+l,mA{F) 

iln+i (m. A) {F) = — — 

AQ,n+l,mA{i-) 

with the collection of integral operators Qn+i,mA from ^[o,n] into ^[o^^+i] defined below 

Qn+iAa, d{a', b)) := 5a{da') v^+^{db) r?f ) (g%) where 9% := G^'\ 
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with B = mA. In the above display d(a',b) = da' x db stands for an infinitesimal neigh- 
borhood of the point {a',b) G ^[o,n+i]) with a = (ao,...,a^) G ^[o,n] and b £ An+i, and 
a = (ao, . . . , On) G -^[o,n]- It is important to point out that 



Bn ■= 7n(l) X An =^ Bn+l = BnQn+1, 



Bn 



Notice that the flow of measures {Bn)n>o satisfies the same type of equation as in 
with the a total mass evolution of the same form as ()1.3p : 

Bn+l{l) = Bnil) An {Gn,Bn) with ^n/mA := j I'n+lidb) Q^^^^^a 

Q„+i,B„(F)(a) = / un+i{db) r/W (g%}^ F{a,b) 



[Qn+l,B{F)-Qn+l,B'{F)] (c 



Un+iidb) 



If we set B = mA and B' = m'A' then condition {H2) is met as long as 

Vi'^ - Vi'^ I < c(n) \m -m'\+ j \[A - A']{^)\ ^%,{d^) 

for some collection of bounded measures ^n^^' ^i^^n) such that f osc((/7) < 



6 {^n^j, for some finite constant 6 {^T^n' j < 00, whose values do dot depend on the pa- 
rameters {m,A) £ {In X V{An))- Under the assumptions (|4.4p . we have 

g%{x) = a^J:HB) G(f)(x) 
for some collection of parameters an\B) satisfying 

a(f)(S) - a(f)(i?')| < c(n) \m -m'\+ j \[A - ^]{^)\ ^%{d^) 



This condition is clearly satisfied for the PHD model discussed in ()4.6p . as long as the 
functions hniun) + dngn{-,yn) are uniformly bounded from above and below. 
For instance, for b = yn ^ {c, c'} we have 



a^\B) 



hr. 



+ j B{da) 4^\dngn{.,b)) 

In this case, we can check that 

- aW(i?')| < c{n)\[B - with (^(f)(a) := r/W(d„5n(., 6)) 

In the same way, we show that the condition (Hi) is also met for the PHD model. This, by 
construction of A^f we find that 



A^ = nn (7^l(l),<-i) + 



N 



N 
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with some local sampling random fields satisfying p.6p . Notice that 
Qn+i (m, A) = ^n„,^A (^) Mn+i,mAia, d{a', b)) 
with the collection of potential functions 

'Hn,mA{a) := Qn+l,myl(l)(«) = vii^ {Qn,mA) 

and the Markov transitions 

KA ( M ' Qn+i,^A{a.d{a',b)) , , . . ^"+1^'^^) ^""'^ (^"%) 
Mn+i,mA{a,d{a ,b)) := — — ^ -— — — = da{da) 



Qn.l,r.A{l){a) ' !.n.l{db'),^:\GZL) 



4.3 Mixed particle association models 

We consider the association mapping 

: (m, A,r/) G (]0, oo[x^[o,n] x ^ J^„+i(m, ^, r/) G P(^[o,„+i]) 

defined for any (m, A) G (]0, oo[x^[q „]) and any mapping 7/ : a G Supp(A) i— )• 7?*^"^ G Pa{En) 

by 

By construction, for any discrete measure A G 'P(^[o.n-i])! and any mapping a G Supp(^) i— 
'qio) g ^{En-i), we have the formula 

[rn,j A{da) r^^'^^)^ = j (m,A,r,^-^) {d{aM (r/(")) 
We also mention that the updating-prediction transformation defined in (|4.5p 

4^^^ (^) = (^) ^'^^ = ^^"J with 4^;) = 5i!)i,^MW (4.8) 

In the above displayed formula 5^"^ stands for some updating Markov transition from En-i 
into itself satisfying the compatibility condition riS^^\ ^ = ^g^ -^ (v)- 

We let ^q'^ = Xli^i '^aj ' be the empirical measure associated with independent and 
identically distributed random variables (aQ)i<i<Ar with common distribution Aq. For any 
a G we let 

,0^ := /<(da) and ry^^'^ = ^ ^ 



• 1 

1 = 1 



with the empirical measure 1]^'^ associated with A^' random variables S^^^ = (^^^''^^ 
with common law rj^f^ . We further assume that 7o(l) is known and set 

7o^:=7o(l)r?o^ and if {!) := 7^ {!) {G,^,n) 

RR n° 7376 



Stability and Approximation of Branching Distribution Flows 



39 



It is readily checked that the fluctuation random fields given below 

,.,(a,Af') /T77 / ia,N') (a)\ 

= ViV' [rj}^ -Vo ) 
satisfies ()1.6p . with N = N' , for any given a G Aq. Using the fact that 



{a,N') 




we conclude that 



with some local sampling random fields satisfying the same estimates as in (|1.6p by 
replacing l/\/iV by the sum (^l/y/N + l/\/]V'^. 
Using ()4.8p . for any ai = (ao,ai) we find that 

We let = Y^i^=i ^a,\ occupation measure associated with N conditionally inde- 

pendent and identically distributed random variables a\ := (oq 1, Oi 1) with common law 

f^i (70^(1), 

In the above displayed formula r?o''^ "* stands for the mapping uq G Aq ^ f?o"°'^ ^ ^ 'P(-E'o)- 
We consider a sequence of conditionally independent random variables with 
distribution /C^'^'J^ ^, (4""'^^' ^^^h 1 < j < iV', and we set 



TV' 

{(ao,ai),Af') _ 1 e „Af ._ /" /jAfrWa^ 

^1 ~ 2^ c'^[(''o.°i)j] and .- / yl^ (daj r^^ 



Arguing as before, for any given ai := (ao,ai) G Supp(^]^), the sequence of random fields 
satisfies (11.61). with N = N'. Thus, we conclude that 



^1 (7o^(l),<,r/i-^'^) (d(«o,ai)) cI>S-) (,(-'^')) + ' Wf 



with some local sampling random fields satisfying the same estimates as in (jl.6p by 
replacing 1/y/N by the sum ^l/\/iV + l/\/]V'^. Iterating this procedure, we define by 
induction a sequence of A^-particle approximation measures 



7;y(l) = 7^-1(1) vtiiG^-i,,^J and r?^ := / (da) r?^^') 
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with the occupation measure = X^i^i ^a.i, associated with N conditionally independent 
and identically distributed random variables aj^ := (og „, a*]^ aj^ „) with common law 

(7,^1(1), A^_i, ^ j. Arguing as above, we find that 



with some local sampling random fields satisfying the same estimates as in (|1.6p by replacing 
l/\/iV by the sum ^l/\/iV + N'^ . As before, the A'^-particle occupation measures A^ 
converge as tends to 00 to the association probability measures A„. 

5 Appendix 

5.1 Proof of corollary 13.61 

For constant mappings s„ = fin+i{l), the mappings ^n^i and 'I'^+i are given by 
'^l+i,uS^) = Sn and $^+i,„„(r/) = ^'^(.) (r^jM^Ji 

with the likelihood function gn "^ and the Markov transitions Mj^^-^ defined in (j2.7|) . Firstly, 

we observe that rn(s„) := sup^. x'eE„ 9n\x) / gn\x'). We also notice that the second compo- 
nent mapping $^+1 does not depends on the parameter m„, and it induces a Feynman- 
Kac semigroup of the same form as the one discussed in section 13.21 

Under the premise of the proposition, the semigroup of associated with the Markov 
transitions M„ satisfies the mixing property stated in the l.h.s. of ()3.1ip for some integer 
m > 1 and some parameter ep{m) g]0, 1]. In this situation, we also have that 

M^pi+^{x, .) > e(^)(m) M^J+„(x', .) 
with some positive parameter 

e^;\m)>epim)/ TT r,{sMl) and r„(5„) := M^^li^ (< r„(l)) 

SnQn + (1 - Sn) 

p<k<p+m ^ ' 

To prove this claim, firstly we observe that Mp^p_,_„(x, .) ^ Mp*p^^(a;, .) and 

p<k<p+m p<k<p+m 
(s) — 

with the semigroup Mp^n associated with the Markov transition 

M^'l~-^{x,.)=ap+iMp+i{x,.) + {l-ak+i)'pk+i with Op+i := ^^^^^ 

Using the geometric representation 

^S"(^'-)= n ak\Mp^n{x,.)+ (l-Ofc) J] ai \ -pkMk,n 

\p<k<n j p<k<n \k<l<n j 
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it can be verified that 

As) 



p<k<p+m 



from which we conclude that 



(^) 



>4'^M M;J+„(x',.) with e^;\m)>epim)/ [] rfc(sfc)rfc(l) 

p<fc<p+m 



We end the proof of the proposition combing the proposition l3.5l with the couple of estimates 
presented in ()3.12p and ()3.13p . This ends the proof of the corollary. ■ 



5.2 Proof of theorem 13.71 

The formulae presented in (|2.6|) can be rewritten in terms of matrix operations as follows 



[7n+l(l) , 1 -7n+l(l)] = [%{!) , 1 -7n(l)] 



/i„+l(l) l-^„+i(l) 



and 



[7„(1) , l-7n(l)] 


Vn{9n) 
1 




[7„(1) , l-7n(l)] 


Vn{9n) 
1 




1 " 
1 



[7n(l) , l-7n(l)] 

With a slight abuse of notation, we set 

i?„ := [7„(1) , 1 - 7n(l)] := [7n(l) , 1 - 7n(l)] and 1 
We also denote by Ain+i,r}n and Vn^rjn the stochastic and the diagonal matrices defined by 



M 



n+l,rin 



and T). 



Vnign) 

1 



(5.1) 



In this notation, the above recursion can be rewritten in a more compact form 



'dn+1 = '&n Mn+l,r^r. and -dr. 



i9. 



n+l 



with the product of matrices Qn+i,r?„ = ^^n,r?„A^n+i,»7„- 

Vn G Ip{C [0, 1]) [^InA^), 1 - ^inA^)] 
with the matrix semigroup 



[u, l-u] Qp^n,u 



[u, 1-u] Qp,n,iy(l) 
Qp,n,u — Qp+l,UpQp+2,Up^i ■ ■ ■ Qn,u„-i 
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These semigroups are again of the same form as the Feynman-Kac models discussed in 
section [3^2] with a two point state space. When G]0, 1[ and < s~ < < 1, we 

have for any n > and any i, i' ,j £ {1, 2} 

Mn+i,uAhj)>enMn+i,uAi',j) and sup g^+i^^^ W(^) < s'^{g) 

The first assertion is a direct consequence of the proposition 13.51 with the couple of estimates 
presented in ()3.12p and (|3.13p . 

Using (j2.3p . we find that induces a Feynman-Kac models of the same form as 

the one discussed in section 13.21 More precisely, we have that 

with the potential functions Gn,m„ and the Markov transitions Mn-\-i,mn defined in ([27 
and ([23]) . Notice that 

sup — — < dn{sg) 



x,x' 



and for any x G En and any n > 

6n{sg)~^ .) < M„+i,^„(x, .) < 5n{sg) M-^^^^^^{x, .) 

with the Markov transitions defined as M^_^-^ by replacing the functions {sn, g-n) 

by their lower bounds {s~,g~). To prove this claim, we use the fact that for any positive 
function / we have 

dMn+l,mAf) ^ rUng^Sn + (1 " mngnSnMn+i (/) + (1 - mn)^n+l (l)7^n+l (/) 

^^n+l,m„(/) "ln5nSn + (1 - Ulngn Sn + (l - nin) fJ,n+li'i^)'Pn+lif) 

and the two series of inequalities 

rUngnSn + (1 " m„)/U„+i(l) 

and _ 

^ ^ mngnSnMn+lif) + (1 - mn)^n+l (Ij/^n+l (/) ^ ^ j^^^^ 

~ m„5t;rs;;: + (1 - m„)^„+i(i)7i„+i(/) 

With a slight abuse of notation, we write Mp^„, and respectively M~j^, the semigroup 

associated with the Markov transitions Mn+i,m„, and resp. M~_)_-|^^^. Using the same 
argument as in the proof of corollary 13.61 it follows that 

M~p+m{x, .) > ep(m) M-p^^{x', .) 

from which we conclude that 

) with ep{m) > ep{m) 5p+k{sg) ^ 

0<fc<m 
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using proposition 13.51 with the couple of estimates presented in (|3.12|) and (|3.13p . we check 
that ()3.3p is satisfied with 

[n / mj — 1 
fc=0 

and some parameters 

0<A;<m 0<k<m 

and 

Pp{m) := ep{m)~^ 6p+k{sg) < ep{m)~^ J| (5p+fc(s5)^ 

0<k<m 0<k<m 

This ends the proof of the first assertion of the theorem. Next, we discuss condition 
(Cont($)). We observe that 



U vignSn) + (1 - ti)/in+l(l) 



After some manipulations 

= uuZ7+i-u) [^..H (^n) 

+ uu{gS+{^~u) uu'[iL)+ll-u) [^9nH(Sn) -^n+l(l)] [y-u'\{gn) 

Recalling that the mapping 9a{x) = ax /{ax + (1 — x)) in increasing on [0, 1] and using the 
fact that 

with the Markov transition 



S^y{x,dx') = 6Adx') + ( 1 - ^ ^g>'){dx') 



9n{x) V an{x) 

we prove 

5n 



V&,„(z.)(.0-^&,„(z.')(.n)| < ^ |(^-^')5n,.'(«n)| (5.2) 

9n 



and for any u £ In = [m^ , m+] 
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This implies that 

'n.+ l — + + I ,'1 +\ - ^ IT- 

m^gn + (1 - mZ) Qn 



H 4— f— — +7 _ _ , 37 Pn-Mn+l|| I — " 

m„filn + (1 - m^) mnQn + (1 - "T-n ) V^n 



< % [(4-«n) + ll'5n-/^n+l(l) 

Using ()2.3p we also find that 



2 / _ mr]{SngnMn+l{f)) + (1 - m) Hn+l{f) 



mr]{sngn) + (1 - 
It is also readily check that 

/i„+i(l) r/(5r„s„) [^'3„^„(r?)M„+i - 7I„+i] (/) {m - m') 



[^n+i,m(r/) (/) = [^^(^^^j + (1 _ m)^„+i(l)] Kr/(s„<7n) + (1 " 

from which we conclude that 

r„+i < sup <^ ^ < 6Jg) sup <^ \ 

ISngn /^n+U-ljJ L Sn /^n+U-lj J 

This ends the proof of the theorem. ■ 



5.3 Proof of proposition 13.31 

The proof of proposition 13.31 is based on the following technical lemma. 

Lemma 5.1 We assume that the regularity conditions {Lip{^)) and {Cont{^)) are satisfied. 
In this situation, for any p < n, u,u' E Ip, r],r]' £ 'P{Ep) and f G Osci{En) and any flow of 
masses and probability measures m = {mn)n>o £ nn>o-^" '^"'^ ^ (^n)n>o G Y\n>o'^('^n) 
we have the following estimates 



p<q<n 

'^ln,m'iv')if)-^ln,Mif)\<<nl \[V ' V'M ^ln,r,' if , d^) + a 



q,n l^-g 



^2 ,^ ^/ 

p<q<n 

with the collection of parameters oL^, i = 1,2, defined in \3. 



Proof: 

We use the decomposition 



+ E [^\,nAK,qA'^'))-^\-l,n,AK,q-lA'^'))\ 



p<q<n 
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and the fact that 



and 

and 



u — u 



$1 



to show that 

p<q<n 

In the same way, we use the decomposition 

+ H [*km(*Lm'(^'))-^^-l,n,m(^L-l,m'(^'))] 



p<q<n 



and the fact that 



g— l,n,m\ p,g— l,m'v'/ // q,n,m \ q—l,q,m I p,q—l,m' \ I /J/ 

$2 ($2 / /NN ^ ^2 /^2 , / /NlN 

q,n,m\ p,q,m,' \i il q,n,m \ q—l,q,m' I p,q—l,m'\l JjJ 



and 



to show that 

q,n,m\^ p,q^m' \ I >> q—l,n,m\ p,q—l,m'\ l >> 

<«^-i.. K-i-K-il 

Using these estimates we conclude that 



mq-i-m i\ 



p<q<n 
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This ends the proof of the lemma. ■ 

Now we come to the proof of proposition 13.31 
Proof of proposition 13. 3t 

We fix a parameter p > 0, and we let (m„)„>p, {m'n)n>p £ Yln>p and {h'n)n>p, and 
{i^n)n>p S Yln>p'^(^n) be defined by the following recursive formulae 

yqyp m'g = ^l^,_^{m'g_^) and z.^ = 

Vg > p mq = <^l„^_^{mg.i) and Vg = ^>^,m^_, (z^g-i) 
with the initial condition for q = p 

{i^P^'^'p) = iv^v') and {mp,mp) = {u,u') 
By construction, we have 

= '^l,q,m'iv') and Vq = $p,g,„(r/) 

as well as 

rn'g = '^l,gy{u') and nig = <^>l^g^^{u) 
In this case, using lemma 15.11 it follows that 

\[Tl,{m\rj')-rl^im,rj)] (/)] 

<ap,„y" \[v-v'K'^)\ ^l,n,^'if,d^)+ al^\Tlg{m',r]') -Tlg{m,r])\ 

p<q<n 

and 

|rp,„(rn',r?') -rp,„(m,r?)| 

<a-l,n\m-m'\+ ^q,nj |[rp,g(m',r/') - r2 g(m,r/)](99)| 

p<q<n 

with the probability measure ^^ g ,^, ,^, = ^l+i^r^ ^{m' ,rj')- 

Combining these two estimates, we arrive at the following inequality 

|[r2^„(m',77')-r2,„(m,77)] (/)] 




p<r<q<n 
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This implies that 



'\+4,nl \[V-V'M ^InAf^d^) 
p<ri<n 

Our next objective is to show that 

|[r2^„(m',V)-rym,r?)] (/)] 

<<Jm-m'|+/3p^,„y" \[ri-r,']{^)\ e'p,n,Af,d^) 

+ ^^i-'-a ••• ^'■'i'" / \i^l,rii^'^v')-'^l,rii^^V)]i^)\ ^l,ri,m' ,r,' (dip) 

p<ri<r2<---<rk<n 

for any k < {n — p) for some Markov transitions 0^ ^ (/, and the parameters 

k-l 

^p,n ~ bp,n + ^p,ri ^ri,r2 • • • bri,n 

1=1 p<ri<...ri<n 

k-l 

Pp,n ~ ^p,n + ^p,ri ^ri,r-2 • • • Wi,n 

1=1 p<ri<...ri<n 

We proceed by induction on the parameter k. Firstly, we observe that the result is satisfied 
for k = 1 with 

We further assume that the result is satisfied at rank k. In this situation, using the fact 
that 

|[r^^,^K,V)-r2,.,(m,r?)] (^)] 

<b'p,r,\m-m'\+al,^J \[r] - r]']i^')\ ^l,^^^,{ip,d^') 

p<ro<ri 

we conclude that 

\[rl^{m',rj')-rl^im,v)] (/)] 

<a';%^\m-m'\ + Plf j |[r?-r/']M| eJ+;^,^,(/, d^^) 
Rin° 7376X K,ri K,r2 ■ ■ ■ K,n J \[^l,roi^' ,v') - 'rl,roi^,V)]i^)\ '^l,ro,m' 
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with 



p<ri<r2<...<ri^<n 

ok+1 = gk ^ 

l^p,n r'p,n ' / j 

p<Ti<r2<---<Tk<n 



X2 ■ ■ ■ "r^:,n 



^p,r\_ bri,r2 ■ ■ ■ ^rfc,n 



and the Markov transition 
d/c+I c\k+l / -P j„\ ak c\k 



e««„,,,(/,dv) = /3*„e«„,,,(/,d^) 

p<ri<r2<...<ri^<n 



'^p,ri ^ri,r2 ■■■ brf.,n i^p,ri,m',ri'^p,ri,r]')i'^f) 



We end the proof of the proposition using the fact that 



|rp,n("i',??') -rp,n("i,??)| < 



p<q<n 



\m — m 



p<(j'<n 



-.2,2 



|[r/-r/'](^')| 



^p,g,m',ri' ®p,q,r}' 



(V) 



This proof of the proposition is now completed. 



5.4 Proof of theorem 13.111 

For any rj £ V{E) and any u, u' G In, we have 



\u — u 



r(l d)+rdh j yn{dy) 



<\u — u'\ 



r(l - rf) + rd/. {jKT3^^ 



This imphes that condition ()3.2p is satisfied with 



9^ 



[h + dm g )]2 

In the same way, for any rj, rj' G V{E) and any u G In, we have 

K+l,viu)-^i+l,r,'iu)=^dhu jynidy) lh+dur,ig{.,ymh+dur,'{g(.,y))] iV - V') {gi- , V)) 

9^ - 9~ 



< rdhm'^ yn 



[h + dm g 



and the probabihty measure 



nlAdv^) oc I yn{dy) g^(^) s (d^) 



[h + dm g 



9+(a)-9 (a) 
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Now, we come to the analysis of the mappings 

^n+i,«(^) r(l - d)u 1JM + J yn{dy) Wu{v,y) ^g(.,y){v)M + ]I 
with the weight functions 

rdur]{g{.,y)) ( h 



h + duri{g{.,y)) \ h + duri{g{.,y)) 
Notice that 

_ rdm-g-{y) . , . rdm+g+{y) 



h + dm g (y) ' h + drn^g~^{y) 

To have a more synthetic formula, we extend the observation state space with two auxiliary 
points ci , C2 and we set 

we extend the likelihood and the weight functions by setting 

gix,ci) = g{x,C2) = 1 

and 

w~{ci) := r(l — d)m~ < Wu{r],ci) := r(l — d)u < w'^{ci) := r(l — d)m'^ 

Wu{V,C2) = W'^{C2) = W^{C2) := 

In this notation, we find that 

with the collection of Markov transitions My defined below 

Vy ^ {02} My = M and M^^ = JI 
Notice that the normalizing constants y^{wu{r], .)) satisfy the following lower bounds 

y^iwuiv,-)) > y'n(.w-)=r{l-d)m'+yn{w-)+fi{l) 

Wc analyze the Lipschitz properties of the mappings ^^+1^^ using the following decom- 
position 

with the signed measures 
and 

A;+i,„(r?,r?') = yc^J^^^^^^^ J ynidy) [wuiv,y) - Wu{v',y)] {^g^.,y){ri')My - $^+i,„(r/')) 
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Arguing as in the proof of theorem 13.71 given in the appendix (see for instance (j5.2p ). one 
checks that 

|A„+i,„(7?,7?')(/)| 

for some collection of Markov transitions S^, from E into itself. It is also readily checked 
that 



hrdm^ 



1 



\h + m dg {y)f 
This clearly implies that condition ()3.3p is satisfied with 



\{r]-r]'){g{.,y))\ 



r{l-d)m+ 



w^g^ 
9' 



hrdm/yn 



- a~ 



{h + m dg )^ 



We analyze the continuity properties of the mappings u i— t- ^n+i ui^) using the following 
decomposition 

This implies that 

- K+i,u'i^)] if) 



< 



This shows that condition (j3.5p is satisfied with 

1 



r(l -d) + hrd 



[h+dm 3 p' 



U — U 



This ends the proof of the theorem 



r(l -d) + hrd % 



(h + dm g )^ 
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